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LONG-TERM  GOALS 

1-  Show  from  observations  whether  the  dynamics  of  the  multi-year  pack  ice  has  a  key  influence  on 
the  location  of  the  following  summer  MIZ. 

2-  Detennine  the  source  regions  of  summer  sea  ice  in  coastal  zones  using  a  back  trajectory  model 
forced  with  satellite-derived  sea-ice  drift  (Maslanik  et  al.  1995,  Emery  et  al.  1997,  Meier  et  al. 
2000,  Tschudi  et  al.  2010) 

3-  Assess  whether  the  source  region  of  sea  ice  melting  in  peripheral  seas  in  the  GCMs  participating 
in  IPCC  AR5  agree  with  observed  source  region  patterns  from  the  satellite-derived  dataset. 

4-  Compare  Lagrangian  ice  trajectories  in  the  model  with  satellite  datasets. 

5-  Repeat  this  comparative  analysis  for  three  critical  climate  horizons:  a  base  period  with  a  dominant 
perennial  sea  ice  cover  and  two  projection  periods  in  the  21st  Century  (2040-2060  and  2080- 
2080). 

OBJECTIVES 

1-  Reduce  uncertainties  in  future  sea  ice  extent  prediction  from  climate  models  participating  in 
CMIP5.  These  include  temporal  evolution  and  geographical  distribution  of  sea  ice  cover  and 
transport  pathways. 

2-  Improve  our  understanding  of  the  strengths  and/or  limitations  of  GCM  predictions  of  future  ice 
edge  positions. 

3-  Identify  limitations  in  the  ability  of  GCMs  to  simulate  source  regions  for  sea  ice  being  advected 
into  marginal  ice  zones  (MIZs)  and  thus  provide  guidance  for  developing  more  reliable  forecasts 
of  the  future  evolution  of  transports  of  material  between  remote  areas  of  the  Arctic. 

4-  Quantify  GCM  biases  in  the  balance  of  thermodynamically  fonned,  in  situ,  ice  production  as 
opposed  to  advection  gain  and  loss. 

5-  Develop  a  seasonal  forecast  model  for  sea  ice  extent  using  satellite-derived  ice  concentrations, 
PIOMAS  ice  thickness,  RGPS  ice  deformation,  Reanalysis  atmospheric  data  and  Lagrangian  sea- 
ice  back  trajectories  to  estimate  thermodynamic  and  dynamic  (advection)  ice  loss. 
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APPROACH 


We  use  a  Lagrangian  trajectory  model  to  track  the  dynamics  of  the  multi-year  and  first  year  ice 
boundary.  To  this  end,  we  are  able  to  identify  the  source  regions  for  sea  ice  that  melts  and  sea  ice  that 
survives  the  summer  melt  season  in  each  of  the  Arctic  peripheral  seas.  The  Lagrangian  Model  is  forced 
with  weekly  mean  satellite-derived  sea-ice  drifts,  and  monthly-  or  daily-mean  sea-ice  drift  from  the 
Coupled  Model  Intercomparison  Project  5  (CMIP5)  GCMs.  The  ice-motion  vectors  come  from  the 
NSIDC’s  CDR  (Climate  Data  Record)  product,  which  is  derived  using  brightness  temperatures  from 
passive  microwave  radiometers  (SSM/I,  SSMIS,  SSMR,  AVHRR,  AMSRE)  using  a  combination  of 
the  bootstrap  and  NASA-Team  algorithm  (Fowler  et  al,  20??).  The  CDR  data  product  also  assimilates 
buoy  data  when  present  within  a  radius  of  a  few  hundred  kilometers.  In  the  summer,  when  the  cross¬ 
correlation  algorithm  between  brightness  temperature  images  has  difficulties  identifying  features  on 
the  ice  surface,  the  data  set  relies  on  both  buoy  data  and  free-drift  estimates,  the  latter  derived  from 
reanalysis  of  sea  level  pressure  data. 

WORK  COMPLETED 

The  Lagrangian  trajectory  model  forced  with  satellite-derived  or  GCMs  sea-ice  drift  data  is  completed. 
The  model  is  also  being  used  by  two  students  involved  in  the  project:  Patricia  de  Repentigny  (Masters 
student)  and  James  Williams  (PhD  student)  working  on  long-tenn  and  seasonal  forecasting  of  Arctic 
sea  ice  respctively. 

The  analysis  of  the  total  thennodynamic  and  dynamic  sea  ice  loss  for  the  instrumental  record  is 
completed.  We  define  thennodynamic  loss  as  sea  ice  loss  resulting  from  basal  melting  via  ocean  heat 
flux,  or  surface  melting  via  radiative/turbulent  losses.  We  define  dynamic  loss  as  summer  sea  ice 
extent  loss  via  sea  ice  export  through  Fram  Strait  (mainly)  or  sea  ice  extent  loss  via  convergence  of  the 
pack  ice  within  the  Artie  Ocean  proper. 

The  quantification  of  dynamical  pre-conditionning  via  advection  of  the  multi-year  pack  ice  is 
complete.  We  define  dynamical  pre-conditioning  as  the  motion  of  the  multi-year  pack  ice  away  from 
coastlines  of  the  peripheral  seas  of  the  Arctic  Ocean.  This  can  be  accomplished  via  multi-year  sea-ice 
export  through  Fram  Strait  (mainly)  or  (since  the  mid-nineties)  through  convegence  of  the  multi-year 
pack  ice  within  the  Arctic  Ocean  proper. 

Analysis  of  all  factors  responsible  for  the  loss  of  summer  sea  ice  (or  the  sea-ice  extent  minimum)  is 
complete.  These  include,  multi-year  ice  advective  flux  away  from  coastlines  in  winter,  Bering  Strait 
heat  inflow  and  basal  ice  loss,  and  surface  melt  during  the  summer  via  radiative/turbulent  fluxes. 

We  have  added  features  to  the  Ice  Tracker  program  originally  developed  by  Chuck  Fowler  at 
University  of  Colorado  (not  publicly  available).  The  new  Ice  Tracker  will  be  hosted  on  the 
http://thepolarhub.org/  web  site  based  at  Columbia  University  that  was  developed  under  a  separate 
PoLAR  initiative  funded  by  the  National  Science  Foundation  (PI:  Pfirman).  The  Ice  Tracker  is  an 
online  program  that  allows  users  to  plot  the  trajectory  of  a  sea  ice  floe  (backward  or  forward  in  time) 
starting  from  any  ice  covered  area  in  the  Arctic  Ocean  within  the  satellite  era  (1979-2013).  The  new 
Ice  Tracker  allows  the  user  to  select  multiple  trajectories,  and  to  overlay  other  fields  (air  temperature, 
ice  age,  ice  thickness,  etc)  on  the  Lagragian  trajectories.  The  data  can  also  be  saved  for  analysis.  This 
tool  is  useful  for  a  wide  variety  of  applications  including  sediment,  nutrient  and  pollutant  transport  by 
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sea  ice,  seasonal  forecasts,  rescue,  emergency  response,  educational  purposes,  etc.  The  site  is  also 
available  on  the  National  Snow  and  Ice  Data  Center  http://extranet.nsidc.org/IceMotion/. 

RESULTS 

The  project  comprises  three  parts:  1-  establishing  the  link,  from  observations,  between  the  dynamic  of 
the  multi-year  ice  pack  in  winter  and  the  following  summer  sea-ice  extent,  with  a  particular  emphasis 
on  the  minimum  sea-ice  extent,  2-  identifying  the  source  regions  for  multi-year  sea  ice  surviving  the 
melt  season  in  each  of  the  Arctic  peripheral  seas  for  today’s  climate  and  for  future  climate  as  simulated 
by  General  Circulation  models,  and  3-  comparing  the  source  region  and  MIZ  dynamics  simulated  by 
GCMs  and  observed  using  the  instrumental  record. 

Results  from  Part  1  show  a  clear  correlation  (r2  =  -0.72)  between  the  state  of  the  previous  winter 
(Dec-Jan-Feb-Mar)  large  scale  atmospheric  circulation  (characterized  here  by  the  Arctic  Oscillation  - 
AO)  and  the  following  summer  minimum  sea  ice  extent  since  1993  when  the  sea  ice  thickness  in  the 
Arctic  appear  to  have  reached  a  critical  value  for  winter  pre-conditioning  to  have  an  impact  on  the 
following  summer  ice  extent  (see  Figure  1).  A  positive  AO  index  also  correlates  with  sea  ice  export 
through  Fram  Strait,  and  winter  surface  air  temperature  and  sea  ice  concentration  anomaly  in  the 
Western  Beaufort,  Chukchi  and  Eurasian  Seas.  The  dynamic  and  thermodynamic  effects  on  the 
location  of  the  MIZ  are  strongly  coupled.  In  fact  southerly  winds  advect  the  (multi-year)  ice  pack  away 
from  coastlines  but  also  advect  warm  air  northward,  limiting  ice  growth  in  the  winter  and  favoring 
melt  in  the  summer.  Our  preliminary  results  using  backward  Lagrangian  trajectories  of  sea  ice  suggest 
that  the  dynamics  (ice  motion)  in  winter  is  at  least  equally  responsible  for  anomalous  sea  ice  extent  and 
position  of  the  MIZ  as  thermodynamic  effect.  Results  also  show  a  strong  correlation  (r  =  0.8)  between 
the  Bering  Strait  heat  inflow  and  sea  ice  concentration  in  the  Chukchi  Sea  region.  In  fact,  it  explains 
30%  of  the  variance  in  the  residual  signal  not  explained  by  the  AO.  These  results  clearly  point  at  the 
importance  of  the  dynamics  of  the  multi-year  pack  ice  in  the  previous  winter  on  the  following  summer 
sea  ice  conditions.  This  is  in  accord  with  results  from  Feltham  et  al.  (2014)  who  showed  a  high 
correlation  between  the  first  year  ice  extent  prior  to  the  melt  season  (May  and  particularly  June)  and 
the  minimum  sea  ice  extent.  Our  analyses  show  that  the  previous  winter  atmospheric  circulation  in 
turns  defines  the  extent  of  the  first  ice  extent  in  May/June.  This  effectively  extends  the  seasonal 
forecast  period  by  3  months  (end  of  March  as  opposed  to  June). 

During  the  backward  Lagrangian  trajectory  analysis,  we  identified  a  significant  positive  bias  between 
the  background  satellite-derived  CDR  dataset  and  buoy  drift  data  (with  buoy  drift  being  faster  by 
several  cm/sec;  see  Figure  2).  This  results  in  a  trend  in  the  mean  Arctic  sea  ice  drift  that  is  much  larger 
than  previously  reported  in  the  literature  using  ice  buoys  alone.  This  trend  is  in  part  the  result  of  the 
increasing  number  of  buoys  assimilated  in  the  satellite  derived  drift  product  as  time  progresses  from 
1979  to  the  present.  We  are  presently  analyzing  the  three  parent  data  sets  used  in  the  construction  of 
the  CDR  satellite-derived  drift  product  in  an  attempt  to  eliminate  this  bias.  The  presence  of  this  bias 
can  potentially  have  important  consequences.  The  CDR  dataset  is  used  to  derive  the  ice-age  data 
product  of  (Fowler  et  al.,  2004).  An  artificial  trend  in  ice  drift  speed,  will  lead  to  an  exaggerated 
temporal  trend  in  the  loss  of  older  ice  in  the  Arctic.  The  CDR  dataset  is  also  used  by  ice-ocean 
modeling  groups  for  instance  to  calibrate/validate  their  ice  thickness  redistribution  scheme  (e.g. 

Hunke,  2014).  An  unforeseen  goal  of  our  work  is  to  now  re-process  the  three  drift  constituents  to  the 
CDR  dataset,  eliminating  this  bias,  in  order  to  resume  the  back  trajectory  work  for  the  instrumental 
record  time  period  on  a  firmer  basis. 


3 


In  preparation  for  the  Lagrangian  back-trajectory  work  using  GCMs,  we  have  run  back  and  forward 
trajectories  using  code  developed  in  the  course  of  this  proposal  and  12-hourly  buoy  data  from  the 
International  Arctic  Buoy  Program.  The  results  show  that  the  use  of  weekly  average  sea-ice  drift  data 
to  calculate  the  back-trajectories  lead  to  relatively  small  error  (150  km  for  a  total  drift  distance  of 
25000  km)  even  for  sea-ice  drifts  of  more  than  2-year.  We  are  now  repeating  the  analysis  using 
monthly  mean  sea  ice  drift  -  typically  saved  by  GCM  modeling  group  participating  in  CMIP5). 
Depending  on  the  error  statistics  using  monthly  mean  drift  data,  we  will  use  the  daily  or  monthly  drift 
speed  from  the  GCM  to  drive  the  Lagrangian  code  and  map  the  regions  for  the  multi-year  ice  surviving 
the  summer  melt  in  each  of  the  Arctic  peripheral  seas  in  todays  climate  and  the  transitional  period  to  a 
summer  ice  free  Arctic  (2040-2060)  and  a  virtually  ice-free  Arctic  (2080-2100). 

IMPACT/APPLICATIONS 

The  CDR  ice  drift  vectors  and  its  companion  data  set,  the  ice  age  data,  are  widely  used  in  the 
community  for  model  validation  and  Arctic  climate  change  assessment.  The  re-processing  of  this 
dataset  will  eliminate  potential  errors  associated  with  a  bias  in  the  satellite-derived  sea  ice  drift  when 
compared  with  the  buoy  data. 

The  Lagrangian  trajectory  model  allow  for  a  clear  attribution  of  summer  sea  ice  loss  between  dynamic 
preconditionning  of  the  system  in  the  winter  and  dynamic  versus  thermodynamic  effect  during  the 
summer. 

The  IceTracker  is  now  publicly  available.  This  new  tool  will  be  publicly  available  online  shortly.  It  can 
be  used  for  education,  environmental,  and  scientific  purposes.  See  Work  Completed  section  for  details. 
Impact:  see  http://blogs.ei.columbia.edu/2013/12/10/explore-the-arctic-ocean-with-icetracker/ 

RELATED  PROJECTS 

See  last  year  for  complete  list. 

Energy  dissipation  in  viscous-plastic  sea-ice  models:  In  viscous-plastic  (VP)  sea-ice  models,  small 
defonnations  are  approximated  by  irreversible  viscous  deformations,  introducing  a  nonphysical  energy 
sink.  As  the  spatial  resolution  and  the  degree  of  numerical  convergence  of  the  models  increase,  linear 
kinematic  features  (LKFs)  are  better  resolved  and  more  states  of  stress  lie  in  the  viscous  regime. 

Energy  dissipation  in  this  nonphysical  viscous  regime  therefore  increases.  We  derive  a  complete 
kinetic  energy  (KE)  balance  for  sea  ice,  including  plastic  and  viscous  energy  sinks  to  study  energy 
dissipation.  The  main  KE  balance  is  between  the  energy  input  by  the  wind  and  the  dissipation  by  the 
water  drag  and  the  internal  stresses  (dissipating  87%  and  13%  of  the  energy  input  on  an  annual 
average).  The  internal  stress  tenn  is  mostly  important  in  winter  when  ice-ice  interactions  are  dominant. 
The  energy  input  that  is  not  dissipated  locally  is  redistributed  laterally  by  the  internal  stresses  into 
regions  of  dissipation  by  small-scale  deformations  (LKFs).  Of  the  13%  dissipated  annually  by  the 
internal  stress  term,  93%  is  dissipated  in  plastic  friction  along  LKFs  (14%  in  ridging,  79%  in  shearing) 
and  7%  is  stored  as  potential  energy  in  ridges.  For  all  time  and  spatial  scales  tested,  the  frictional 
viscous  dissipation  is  negligible  in  the  KE  balance.  This  conclusion  remains  valid  regardless  of  the 
degree  of  numerical  convergence  of  the  simulations.  Overall,  the  results  confirm  the  applicability,  from 
an  energetical  point  of  view,  of  the  VP  approximation. 
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